Disorder-Induced Anomalous Kinetics in the A + A — > Reaction 
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We address the two-dimensional bimolecular annihilation reaction A + A — > in the presence 
of random impurities. Impurities with sufficiently long-ranged interaction energies are known to 
lead to anomalous diffusion, (r 2 (t)) ~ £ 1-<s , in the absence of reaction. Applying renormalization 
group theory to a field theoretic description of this reaction, we find that this disorder also leads to 
anomalous kinetics in the long time limit: c(t) ~ t s_1 . This kinetics results because the disorder 
forces the system into the (sub) diffusion controlled regime, in which the kinetics must become 
anomalous. 
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Surface reactions show a variety of complex spatial 
and temporal patterns. Simple systems, such as oxi- 
dation of CO on single crystal Pt(llO), show surpris- 
ingly rich behavior ffl , ranging from spirals and standing 
waves to chemical turbulence [||J3) . Such behavior results 
because two dimensions is the upper critical dimension 
for many surface reactions, and so collective fluctuations 
control the dynamics. Toussaint and Wilczek, and Kang 
and Redner first noticed the diffusion-controlled kinetics 
that can arise from reactant microphase separation ^-|) . 
Field theoretic techniques were developed that rigorously 
showed two dimensions to be the upper critical dimen- 
sion for the annihilation reactions A + A — » [[7]j|] and 
A + B — > || . Adsorption effects have been addressed 
within the field theoretic framework [^0|. Real systems, 
of course, possess many defects, and randomness in ad- 
sorption energies can lead to a variety of phases observed 
at steady state [ O . 

Defects can dramatically affect the diffusion of parti- 
cles if correlations in the potential field are sufficiently 
long ranged. For a single ion diffusing in disorder that is 
charged, the appropriate form of the correlation function 
at long-wavelengths is fl2|| 



Xw(k) = l/k 2 



(1) 



For this model, the upper critical dimension is two. In 
two dimensions and below, subdiffusion occurs [p^|Jl4[ . In 
this anomalous regime the mean-square displacement in- 
creases sublinearly with time, (r 2 (t)) ~ t x ~ s . The diffu- 
sion exponent, 1 — 5, depends on the strength of disorder; 
it can be found exactly [|l5|-^lf . 

These same defects can have an effect on reacting 
species. We consider species that diffuse in this quenched, 
random potential and which also react according to the 
law A + A — > 0. The diffusing species, however, do not 
interact in any other way. 

The type of physical system that we have in mind is a 
reaction that occurs on the surface of a crystalline ionic 
lattice. The substrate lattice has dislocation line pairs, 
which form line vacancies or line interstitials. These de- 
fects are immobile and generate a random, quenched elec- 
trostatic potential on the surface, represented by Eq. (j^). 



On the surface, the reaction A + A — > occurs, where A 
is a reactive ion. The 1 jr interaction between the ions is 
technically irrelevant, and so it can be ignored. So as to 
maintain surface charge neutrality, however, we may wish 
to include additional counter ions, which can desorb. 

Alternatively, our results can be viewed as approxi- 
mately describing the annealing of two-dimensional topo- 
logical defects, such as line dislocations in solids. Topo- 
logical defects have logarithmic interactions, due to a 
slowly decaying elastic strain field (see |E| for a re- 
view). If some of the defects were pinned and unreactive, 
the remaining defects would diffuse and combine in the 
quenched, random potential described by Eq. ([!]). The 
logarithmic interaction between the topological defects 
prevents microphase separation, and so the kinetics of 
this two-species reaction is similar to that of the A + A 
reaction Q. 

In this Rapid Communication we analyze a field theo- 
retic formulation of the two-dimensional annihilation re- 
action A + A — > in the presence of quenched disorder. 
The concentration of A is initially Poissonian, with aver- 
age density no- No creation or unimolecular annihilation 
is allowed. We search for the asymptotic decay law at 
long times using renormalization group (RG) theory. A 
field theory is derived by identifying a master equation, 
writing the master equation in terms of creation and an- 
nihilation operators, and using the coherent state rep- 
resentation [0J^]. The random potential is incorporated 
with the replica trick jl4]]. The concentration of A at 
time t, c(x, t), is given by 



c(x, i) = Urn (<z(x, t)) , 



(2) 



where the average is taken with respect to exp(— S), with 
the action S = S + Si + S 2 + S 3 , 

S = [ d d x [ f dta a (x,t) [d t - DV 2 + 6(t)] a a {x,t) 



Si 



,/''x / f dt 





Aia a (x,t)a 2 1 (x,t) 



+A 2 a 2 (x,t)a 2 (x, t) 



1 



S 2 = -no J d d xa a (x,t) 
(3 2 D 2 



S, = 



dtidt^ 



kik 2 k 3 k4 



x(27r) d 5(k 1 +k 2 + k3 + k 4 ) 

xa ai (ki, ti)a 01 (k 2 ,ti)a a2 (k 3 , t 2 )d a2 (k 4 , t 2 ) 

xkx • (ki + k 2 )k 3 • (ki + k 2 )x™(|k 1 + k 2 |) . (3) 

Summation is implied over replica indices. The notation 
J k stands for J d d \t/(2-K) d . The upper time limit in the 
action is arbitrary as long as tf > t. The term 5*0 rep- 
resents simple diffusion, without an external potential. 
The delta function, often left out by convention, enforces 
the initial condition on the free field propagator. The 
term Si comes from the reaction terms. The parame- 
ters are related to the conventional reaction rate, k, by 
Ai = 2A 2 = k. The term S 2 comes from the random, 
Poissonian initial condition. The term S3 comes from 
averaging the concentration over the random potential. 
The potential is assumed to be Gaussian, with zero mean 
and correlation function given by Eq. (Q). The inverse 
temperature is given by (5 — 1/(&bT). 

An approximate solution to Eq. (||) can be derived by a 
saddle point approximation to the action (||). The result, 
before the average over the random potential is taken, is 

d t c v = DV 2 c v + [3DV ■ (c v Vv) - Xicl 
c.„(x, 0) = n . (4) 

For constant potential, this equation has the solution 
c„(x,i) = l/(Ait+l/no). 

An alternative, exact expression for the concentration 
can be derived by performing a Hubbard-Stratonovich 
transformation on (^|) and integrating out the field a: 

d t c nv = DV 2 c nv + (3DV ■ (crfyVv) - Xic 2 



c vv (yL, 0) = n , 



(5) 



where the real, Gaussian, random field 77 has zero mean 
and variance (r?(x,t)r?(x',t')> = 2X 2 5(-K-x')S(t-t'). The 
physical concentration is given by averaging the solution 
over the random field 77. From this representation, we see 
that mean field theory is exact when A 2 = 0. It is for this 
reason that diagrams to all orders in Ai can be summed 
in this limit (see, for example, ||). 

We apply renormalization group theory to the action 
(|J) to take into account the effects of nonzero A 2 . Time 
ordering prevents a term of the form da from being gen- 
erated, and all other relevant terms are already present 
in Eq. (|J). The limit N — > turns out not to mat- 
ter, as time ordering cancels the same terms that this 
limit does. We follow the flow equations until A 2 is small 
enough so that we have entered the perturbative regime. 
In this regime, we can match the flow equations with 



the solution to Eq. (Q). We integrate over momenta in 
the range A/6 < k < A and then rescale the fields by 
a (bk,b- z t) = o(k,t)/a and a'(6k,6" 2 i) = a(k,t)/a. To 
achieve a fixed point, and to keep the time derivative in 
So constant, we set a = 1, a — b d . We determine the 
dynamical exponent by requiring that the diffusion co- 
efficient remain unchanged. The flow equations in two 
dimensions, to one loop order, then become 



din no 


= 2 




dl 






din Ai 


A 2 




dl 


2ttD 


4-7T 


din A 2 


A 2 




dl 


2ttD 


4-7T 


dP 2 ! 
dl 


= . 





The dynamical exponent is given by 

/3 2 7 



z = 2 



4tt 



(6) 



(7) 



These flow equations reproduce known anomalous scal- 
ing in the cases of no reaction or no disorder. They are 
integrated to a time such that 



t(F) = te 



. - [' z(l)dl 



tn 



(8) 



The matching time, to, is chosen to be on the order of 
47r 2 / (A 2 D) so as to be within the range of validity of both 
RG scaling and mean field theory. Mean field theory is a 
good approximation at this time because A 2 (7*) is small. 
(An expansion in A 2 generates an expansion in to that 
leads to subdominant scaling.) For short times, particles 
see only the local value of the potential, and the effective 
diffusivity is given by the bare value. (An expansion in /3 
generates an expansion in to that leads to subdominant 
scaling.) We can, therefore, ignore the potential term in 
Eq. (|J) for short times. Matching with mean field theory, 
we find the mean square displacement is given by 



(r 2 (i(P),0) = 4D*(Z*) 



(9) 



in the absence of reaction, and the concentration profile 
is given by 



1 



c(i(Z*),n 1/nQ (i*) + Xl (l*)t(l*) 



(10) 



in the absence of random forces. The observed values are 
related by scaling: 



(r 2 (t)) =e 2l (r 2 (t(l),l)) 
c(t) = e- 2l c(t(l),l) . 



(11) 



With no reaction, the one loop flow equations prove 
to be exact Jl5|-^l[. The disorder strength is not renor- 
malized to any order. The dynamical exponent is greater 



2 



than two, and so we have subdiffusive behavior. Match- 
ing at I = I*, we find the mean square displacement at 
long times is given by 



(r 2 (t)) ~ ADt 



with 



5 = 



1 



/3 2 7. 



(12) 



(13) 



We have a continuously variable exponent that depends 
on the strength of disorder. The prefactor is nonuniver- 
sal, since the matching time to explicitly enters. 

With no external potential, the flow equations become 
asymptotically exact at long times [Mil . The rate con- 
stants decay slowly to zero, with the ratio A2 (Z) /Ai (Z ) re- 
maining constant. The dynamical exponent is two, since 
there is no field renormalization. Matching with mean 
field theory at I = /*, we find for long times 



c(t) 



Inft/to) 



(14) 



We see the characteristic universal logarithmic correction 
to the mean-field kinetic equations. There is a subdomi- 
nant, nonuniversal 1/t decay. 

The presense of disorder dramatically affects the reac- 
tion. The rate constants no longer decay to zero, but to a 
nonzero fixed point. The ratio A2 (0 / Ai (Z) again remains 
constant. To one loop, we find the reaction terms do not 
renormalize the diffusivity, disorder strength, or dynami- 
cal exponent. Matching with mean field theory at I = I*, 
we find 



c(t) 



— (- 

with the fixed point given from Eq. as 
XI = 3Df3 2 -f . 



(15) 



(16) 



We have, therefore, a nonuniversal decay exponent that 
depends, again, on the strength of disorder. The expo- 
nent is precisely the inverse of that for the pure diffusion 
case, Eq. (|l2|). The prefactor is also nonuniversal, since 
the matching time to again enters. 

We see that to one loop, the reaction terms leave 
the diffusion process unaffected. Physically, the disor- 
der traps diffusing species in wells of deep potential en- 
ergy. Normal diffusion between these wells leads to the 
observed, overall diffusivity. Progressively deeper wells 
are encountered at longer times, and so anomalous dif- 
fusion occurs. For diffusion limited reactions, the effec- 
tive reaction rate is proportional to the effective diffusion 
constant. Since we have D e ff(i) ~ D(t/to)~ d , we should 
expect k e ff oc (t/to)~ s ■ This is exactly the behavior ob- 
served in Eq. (|i1f). 



When two reactants are attracted to the same well, 
reaction almost always occurs, since the escape time is 
much greater than the reaction time. In other words, the 
reaction becomes (sub) diffusion limited, and the bare re- 
action rate should not enter the decay law. Three species 
entering a well is much less probable at low densities than 
two, and so this case should be irrelevant at long times. 
Reactants can encounter each other outside the wells. 
Reaction outside the wells would show up as renormal- 
ization of the disorder strength and dynamical exponent 
in a two-loop calculation. We have checked that these 
terms are not generated to two-loop order. So the dy- 
namical exponent is unchanged by the reaction, and the 
bare reaction rate is not present in the decay law. From 
this argument, we expect the dynamical exponent and 
the flow equation for the disorder to remain unchanged 
in higher order calculations. 

The diffusion terms, on the other hand, do affect the 
reaction process. Since the reaction becomes diffusion 
limited, the diffusion terms actually control the reaction 
process. For this reason, we find a nonzero fixed point 
for the effective reaction rate. The effective reaction rate 
should be a function of the strength of disorder, since 
it is related to the rate at which species diffuse between 
wells. The rate (at I = I*) is proportional to the effective 
diffusion coefficient divided by a characteristic distance 
squared. This is exactly the behavior that we observe in 



Eq. (16), with /3 2 7 an inverse distance squared. Higher 
order calculations may, of course, lead to modification of 
the numerical value of the fixed point for large /3 2 7- 
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